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ABSTRACT 

Constraining the spatial and thermal structure of the gaseous component of 
circumstellar disks is crucial to understand star and planet formation. Mod- 
els predict that the [Ne II] line at 12.81 /im detected in young stellar objects 
with Spitzer traces disk gas and its response to high energy radiation, but such 
[Ne II] emission may also originate in shocks within powerful outflows. To dis- 
tinguish between these potential origins for mid-infrared [Ne II] emission and to 
constrain disk models, we observed 32 young stellar objects using the high reso- 
lution (R~30000) mid-infrared spectrograph VISIR at the VLT. We detected the 
12.81 /im [Ne II] line in 12 objects, tripling the number of detections of this line 
in young stellar objects with high spatial and spectral resolution spectrographs. 
We obtain the following main results: a) In Class I objects the [Ne II] emission 
observed from Spitzer is mainly due to gas at a distance of more than 20-40 AU 
from the star, where neon is, most likely, ionized by shocks due to protostel- 
lar outflows, b) In transition and pre-transition disks, most of the emission is 
confined to the inner disk, within 20-40 AU from the central star, c) Detailed 
analysis of line profiles indicates that, in transition and pre-transition disks, the 
line is slightly blue-shifted (2-12 km s" 1 ) with respect to the stellar velocity, and 
the line width is directly correlated with the disk inclination, as expected if the 
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emission is due to a disk wind, d) Models of EUV/X-ray irradiated disks repro- 
duce well the observed relation between the line width and the disk inclination, 
but underestimate the blue-shift of the line. 

Subject headings: stars:formation, infrared: stars, protoplanetary disks, circum- 
stellar matter, line:profiles 



Introduction 



Disks orbiting young stars (1-10 Myr) are sources of accreting material and sites of 
nascent planetary systems. Observations that can establish the physical conditions and evo- 
lution of the gaseous component within young circumstellar disks are essential if we are to 
understand the magnetospheric accretion processes that determines a star's zero age main 
sequence mass, protoplanetary disk structure and evolution, and the processes involved in 
planet formation. Despite their relevance to all of these important open issues, the physical 
and chemical properties of gaseous circumstellar disks are still poorly understood. Submil- 
limeter observations of molecular lines provi de insight into the chem i stry and the physics 



of th e outer disk regions (r>100 AU) (e.g. iDutrey et al.l Il994 Il997l ; Kastner et al.l Il997 



20081 ). but the properties of the gas within 20-30 AU have been only recently investigated 
thanks to the mid- infrared IRS spectrograph on board the Spitzer Space Telescope, e.g., by 
the discovery of emission lines from se veral species including organic molecules and water 
f lCarr fc Naiitall2008[ ISaTvk et al.ll2008h . 



One of the most interesting mid-infrared spectral features is the [Ne II] line at ^.Sl - 



um (IPascucci et al.l 120071 ; lLahuis et al.l l2007t lEspaillat et al.l l2007bl ; iFlaccomio et al.l 12009 



Giidel et al.ll2010l ). The importance of the [Ne II] line is twofold: it traces warm gas (T~5000 
K) and, due to the high first ionization potential of neon (21.56 eV), it is a good proxy for 
the influence of stellar extreme ultraviolet (EUV) and X-ray radiation on the disk. The 
presence of strong [Ne II ] emiss ion at 12.81 /j,m from protoplanetary disks was first pre- 
dicted by iGlassgold et al.l ( 120071 ). who suggested that neon is ionized by K-shell absorption 
of X-ray photons at en ergy E>0.9 keV. Thi s hypothesis has been refined by sophisticate d 
irradiated-disk models (IMeijerink et al.ll2008l : lErcolano et al.ll2008l : lErcolano &: Owenll2010l ). 
which predict a correlation between the stellar X-ray luminosity and the total [Ne II] line 



1 Based on observations collected at the European Organisation for Astronomical Research in the Southern 
Hemisphere, Chile, with the mid-infrared spectrograph and imager VISIR (Prog. ID 083.C-0883, 084.C-0211 
and 085.C-0860) 
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luminosity. It has also been suggested that X-rays may determine disk dissipation by heat- 
ing the circumstellar gas to t emperatures 10 3 — 10 4 K and trigger ing its photoevap- 



oration (lErcolano et al.ll2008l ; lErcolano &: Owenll2010l ; lOwen et al.ll2010l ). Neon ionization 



and dis k photoevaporation could also be induced by EUV photons produced by accretion 
shocks (IHollenbach fc Gortill2009l ; lAlexanderl 120081 ) . However, EUV photons are more easily 
absorbed by neutral hydrogen, so they are less likely than X-rays to reach and affect the 
circumstellar disk, esp ecially early in the C l ass II phase when inflow and outflow from the 
inner disk are densest ( lErcolano et al.ll2009l ; lOwen et al.ll2010l ). 



A detailed comparison between models and observations is required to clarify this com- 
plex theoretical scenario. Several surveys of the [Ne II] emission from young stellar objects 
(YSOs) have been carried out in the last 5 years, u s ing low spectral and spatial resolution data 
obtained with Spitzer /IRS dPascucci et al. 2007 ; Espaillat et al. 2007b ; Lahuis et al. 2007 



Flaccomio et al. 2009; Giidel et al. 



2010h. A we ak correlation between X-ray and [Ne II] 



luminosity has been found by iGiidel et al.l (120101 ) , who investigated a sample consisting of 
92 pre-main sequence stars. Furthermore, [Ne II] emission was found to be stronger in Class 



I YSOs and in sources d riving protostellar jets than in more evolved YSOs (jFlaccomio et al. 



20091 ; IGiidel et al.l 120101 ) . The strong [Ne II] emission observed in embedded YSOs driving 



jets - orders of magnitude larger than predicted by models of emission from disks - is un- 
likely produced by irradiated gas within the inner disk, but it may be produced eithe r by 
shocks, in the circumstellar envelope or in protostellar jets (IHollenbach fc Gortil 120091 ). or 
i n the launching re gion of magnetically driven outflows irradiated by stellar X-ray emission 
Jshang et alJboioh . 



The Spitzer /IRS data, with their limited spectral and spatial resolution, provide limited 
information on the velocity and the spatial structure of the emitting gas. High spectral and 
spatial resolution observations of the [Ne II] line from the ground can be used to decouple 
disk and jet emission. For example, using VLT/VISIR, clear evidence o f [Ne II] emission 



from protostellar jets has been discovered in the triple system T Tau (Ivan Boekel et al. 



20091 ). The observations show that the predominant [Ne II] emission component is shifted to 
a velocity compatible with the jet motion and is extended (~ 1.1 ) along the jet direction. 
However, a disk origin for the [Ne II] emission has been suggested by hi gh spectral resolutio n 
observations of the ne arby classical T Ta uri stars (cTTSs) TW Hya (IHerczeg et al.l 120071 ). 
AA Tau and GM Tau (INajita et al.ll2009f ). In all of these stars the line appeared broadened 
(Av = 21 — 70 km and centered near the stellar velocity. The absence of detectable 
blue-shift suggested that the emission is produced from a static disk atmosphere, but given 
the low signal-to-noise ratio of this data, it was not possible to rule out the presence of 
emission from a slow (v = 1 — 10 km photoevaporative wind. 
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Evidence of [Ne II] emis sion arising from both a pr otostellar jet and a photoevaporative 
disk wind has been found by lPascucci Sterzikl (120091 ). who observed 4 optically thick disk 
objects (Sz 73, Sz 102, HD 3700A, VW Cha) and 3 transition disk objects (TW Hydrae, 
T Cha and CS Cha). Transition disk objects have no excess in the near infrared and a 
strong excess in the far infrared, which is inter preted as a signature o f a few AU central 
gap w ithin the disk (see I Williams fc Ciezall201ll and references therein). iPascucci fc Sterzik 
(120091 ) detected [Ne II] emission from all the transition disk objects, finding the line slightly 
blue-shifted (Av ~ 4 — 6 km s _1 ) and broadened, with a line width correlated with incli- 
nation angle. These characteristics appear consistent with model profiles expected from a 



photoevaporative disk w ind triggered by EU V or X-ray stellar emission ([Alexander! 12008 



Ercolano &: Owen! |2010| ). but a larger number of analogous observations are required to 
confirm their results. They detected the emission in only one of the optically thick disk 
objects, Sz 73. Since the emission is strongly blue-shifted (Av = 100 km s" 1 ) and broadened 
(Av = 60 km s _1 ), an origin from fast outflowing material is more likely than a disk origin. 
Additional evidence for a photoevaporative wind comes from the asymmetric [Ne II] line 
profile detected toward the transition disk around TW Hydrae (IPascucci et al.ll201lh . 

With the aim of understanding the origin of the [Ne II] emission from stars at various 
evolutionary stages and with different circumstellar characteristics and to test X-ray and 
EUV irradiated disk models with observations, we performed a high resolution spectroscopic 
survey of the [Ne II] emission from 32 T Tauri stars and YSOs spanning a range of spectral 
energy distributions (SED), with the VISIR spectrograph at the VLT. Here, we present 
the results of this survey, including a comparison with Spitzer /IRS spectra, some of which 
are published here for the first time. In Sect. [21 we describe the observed sample, the 
observations and the data reduction reduction procedure; in Sect. [3] we report the results of 
analysis of the VISIR and the Spitzer data; in Sect. 4, we discuss our results and compare 
them with the predictions of irradiated disk models; and in Sect. 5, we draw our conclusions. 



2. Observations and data reduction 

We observed 32 T Tauri stars and YSOs (2 resolved binaries IRS 43 and CrA IRS5 are 
counted as 4), belonging to several star forming regions and young stellar associations. In 
subsection I2.1[ we describe the properties of the observed stars and the procedures used to 
observe them; in subsection 12. 2} we describe the data reduction process. 
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2.1. Target sample and observations 

All the stars in our sample, with the exception of IRAS 08267-3336, are closer than 
150 pc. We include stars from both very young star forming regions such as p Oph and 
Corona Australis (1-3 Myr), which are characterized by a large number of deeply embedded 
protostellar objects, and older stellar associations like e Cha and /3 Pictoris (6-12 Myr), that 
are mainly composed of apparently diskless (Class III) stars. We consider different regions 
and associations to include in our sample YSOs of different infrared Classes. Furthermore, 
we gave higher priority to YSOs with Spitzer detections of the [Ne II] emission at 12.81 pm. 
T wo objects (T Cha and Sz 73) - whose [Ne II] line was already detected with VLT/VISIR 



by iPascucci fc Sterziki (120091 ) - have been observed again to investigate the variability of the 



emission. 



The full list of the observed targets is reported in Table [fl where we list object coor- 



dinates from the 2MASS catalog (ISkrutskie et al.ll2006l ). parent region, distance (based on 



the parent region and data reported in the literature), heliocentric radial velocities, infrared 
Class, disk inclination, X-ray luminosity in the band 0.3-10.0 keV@ and mass accretion rate. 
All data in Table [1] are retrieved from the literature. Radial velocities for some of the stars 
in the p Oph and Corona Australis star forming regions were not available in the litera- 
ture, so these radial velocities have been estimated from the velocity of the parent cloud 
or nearby stars (see Tabled] for further details). Disk inclinations reported in Table [T] are 
derived by different techniques. Specifically, in 8 YSOs, disk inclinations are derived from 
interferometric submillimeter observations or scattered-light imaging of the disk (with esti- 
mated uncertainties 5-10 degrees). For RU Lupi we assumed the disk inclination to be equal 
to the inclination of the stellar rotation axes (error ~10 degrees), and for T Cha and RX 
J1615. 3-3255 we report the inclination derived from a fit of the spectral energy distribution 
(SED). Errors on the inclination derived from the SED are probably much higher. However, 
in the case of T Cha (inclination angle 75 degrees) other observational signatures - e.g. a 
highly variable optical extinction - demonstrate the presence of a circumstellar disk seen 
nearly edge-on. For all the remaining targets, no information on the disk inclinations are 
available in the literature. 

The infrared Class of our targets, reported in Table[U was derived from analyses reported 
in the literature (references are reported in the last column of Table [T]) of the SED in the 
near-infrared and mid-infrared bands and, in three cases (R CrA-7B, CrA IRS 5 and T 
CrA), of mid- infrared color- magnitude diagrams. Specifically, we classified, as Class I and 



2 Some published Lx values are for slightly different bands and were corrected to account for the band 
differences. 
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Class II, sources with SED slopes in the ranges a > and -2 < a <0, respectively. To 
first approximation, the slope of the SED indicates the presence of a massive protostellar 
envelope (in Class I) or an optically thick disk (in Class II). However , recent detailed studies 



of the dust emission in the mid-infrared and submillimeter band (see Williams fc Cieza 2011 



for a review) suggest that within the Class II are included stars with an optically thick 
disk and stars harbouri ng a more evolved disk, c haracterized by the presence of an inner 
gap. As pointed out by iPascucci fc Sterzikl ( 120091 ). the origin of the [Ne II] emission in the 
latter group of stars is much different than in the former, therefore, following the definitions 
mostly used in the literature, we define two more Classes of objects: the transition disks 
and the pre-transition disks. Specifically, circumstellar disks with signatures of a large (few 
AU) inner hole - i.e. no excess in the near infrared and large excess in the mid infrared - 
are classified as transition objects (Tr), while circumstellar disks, with a SED indicating the 
presence of a gap, which separates the inner d isk from the outer disk are classified as pre- 
transition disks (pTr, see lEspaillat et al.ll2007al for a detailed discussion about the difference 
between transition and pre-transition disks). 

On the basis of these criteria, our sample is composed of 9 Class I, 13 Class II, 6 
transition disks, and 4 pre-transition disks. However, the information that we found in the 
literature for our sample were not homogeneous so, we cannot exclude that with a more 
complete and homogeneous set of data some of the Class II objects may be classified as 
transition or pre-transition disks. 

It is worth noting that a large gap in the dusty component of the disk does not necessarily 
imply the presence of a gap in its gaseous component, where the [Ne II] emission is produced. 
For example, in our sample T Cha, V4046 Sgr and Hen 3-600 are transition disk objects but 



they are also actively acc reting, as strong optical emission lines demonstrate ( jSchisano et al. 



20091 ; ICurran et al.ll20f it ). The presence of these accretion signatures indicates that the hole 
within the dusty component of the disk is filled with gas. 

VISIR observations have been ca rried out during; thr ee different ESO observing periods, 
P83, P84 and P85. We used VISIR (ILagage et al.l l2004h in high resolution, long slit, spec- 
tros copic mode. The ech elle grating is centered on the [Ne II] fine structure line at 12.81355 
fj,m ( lYamada et al.lll985l ) and covers the spectral range between 12.79 and 12.83 fim. We set 

the slit width at 0.4 (~40 AU at a distance of ~100 pc), which coincides with the spatial 

" I 

resolu tion of the instrument when the visual seeing is better than 0.6-0.7 (ISmette fc Vanzi 
20071 ). and a spectral resolution R~30,000 (~ 10 km s -1 in velocity scale). The spectral and 
spatial resolutions of our observations are, therefore, 50 and 10 times higher than [Ne II] 
observations carried out using IRS spectrograph on board Spitzer, respectively (R~600 and 
slit width ~ 4.7 ). For flux calibration and telluric correction, we observed a standard star 
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immediately before and/or after the observation of each target at airmass as close as possible 
to that of the target. The observation of stars requiring long exposure times have been split 
in more than one segment, alternating an observation of the target with an observation of 
the standard star. Table [2] reports date, time, airmass and exposure time of all the observed 
stars and their associated standard stars. During the first run, we also observed the aster- 
oid Kalliope to obtain a pure sky absorption spectru m, and Titan to verify the instrument 
resolution. Our Titan spectrum has been analyzed by lPascucci et al.l ( 120111 ). who measured 
the FWHMs and the centroid position of unresolved lines of C 2 H 2 and C 2 H 6 . The measured 
FWHMs span a very narrow range with a mean value of 10 km s" 1 , as expected considering 
the instrument resolution, while the uncertainties on the line centroids is ~1 km s _1 . How- 



ever, by tests on a set of standard stars and from TW Hydrae data, iPascucci et al.l (1201 if ) 
found that the absolute accuracy of the line centroid is 2 km s _1 due to uncertainties in 
placing the source in the center of the slit. 



We also compiled fl uxes in the [Ne II] lines measured from Spitzer /IRS ( 1 Werner et al. 



2004 ; iHouck et al.l 120041 ) from the literature, and when these were not available, we down- 
loaded, reduced and analyzed data from the Spitzer archive. Spitzer [Ne II] fluxes are 
reported in Table [3J together with the results obtained from the VISIR observations, for a 
direct comparison. 



2.2. Data reduction 



We reduced all VLT/VISIR data using the VISIR pipeline Vers. 3.2.2 (ILundinl 120081 ) 
in conjunction with our IDL scripts. VISIR data are a collection of data cubes, each corre- 
sponding to a nod position. Every data cube contains 2n+l planes, where n is the number 
of the chopping cycles. Odd planes contain the images from the on source chopping cycle 
(Ai), while the even planes contain the average of the differences between the current and 
all the previous on source (Ai) and off source (BA chopping cycle images. Therefore, the 
last plane contains the average of all the Ai — Bi images. 

The first part of the data reduction is executed by the pipeline. Specifically, all the data 
cubes are coupled in pairs, composed of two complementary nod positions. The last planes 
of each pair are summed and divided by the total integration time. At this point, because we 
used standard observing mode, in which the chopping and nodding directions and throws are 
the same, each nodded image contains a double positive beam in the center and a negative 
beam on each side. To form the final science frame, all the nodded images are added, after 
bad pixel cleaning and correction of the image distortion. In a small number of cases, we 
excluded from the final sum some pairs of nod images taken when the background emission 
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was particularly strong and variable, e.g. due to the presence of thin clouds. The standard 
deviation of the noise in the resulting image is calculated from the dispersion of pixel values, 
after the rejection of outliers with an iterative sigma-clipping. Pixels with values less than 
3 times the standard deviation of the noise are classified as noisy. 



The spectrum is extracted by the pipeline using the method developed by iHornd (119861 ). 
Specifically, a weight map is generated by collapsing the spectrum along the dispersion 
direction, normalizing the absolute flux of the one-dimensional image to 1 and setting to 
zero noisy pixels. The resulting profile is then replicated along the dispersion direction to 
cover the full frame. The final spectrum is obtained by multiplying the science frame for 
the weight map and summing along the slit direction. Pixels classified as noisy are used to 
calculate errors (see Fig. [Q. 

After several tests, we realized that this method for spectrum extraction fails in two 
situations: a) when the signal-to-noise ratio is very low and the residual noise is comparable 
to the signal from the star, because the weight map, and therefore the extracted spectrum, 
is significantly contaminated by noise (see the the middle panels of Fig. [T]); b) when we 
observe a spatially resolved binary with the two components aligned along the slit, because 
by multiplying the science frame for the weight map and summing along the slit direction 
we obtain an integrated spectrum of both components, losing information on the emission 
of each component. To avoid these problems, we extracted the spectra of faint sources and 
resolved binary systems by defining our own weight map. The weight map for faint sources is 
the sum of three gaussians, one centered on the peaks of the positive beam and two centered 
on the negative beams, with half and opposite amplitudes with respect to the first one 
(see Fig. [TJ). To extract separately the spectra of each component of the spatially resolved 
binaries, we defined different maps for each component, using the same method described 
above (see bottom panels of Fig. [Q To check this procedure, we applied this method to the 
bright stars of our sample; obtaining the same results as with the standard procedure (see 
top panels of Fig. [p. 

The wavelength calibration is done by cross- correlating an infrared background spectrum 
with a synthetic model spectrum of the atmosphere. The infrared background spectrum 
is extracted from the input files before correcting for chopping and nodding. The whole 
extraction procedure described above is performed in the same way for both the science 
target and the associated standard stars. Finally, the flux calibrated spectrum is obtained 
by multiplying the observed spectrum of the star by a spectral response function. The 
response function is calculated by dividing the observed spectrum of a standard star by the 
flux model (at very low resolution and not including photospheric absorption lines) retrieved 
from the VISIR standard stars catalog. Photospheric absorption lines in the spectra of the 
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standard stars have been eliminated by gaussian fitting. This procedure allow us to correct 
both for sky transmission, which depends on the wavelength - because of telluric lines - and 
on the airmass, and for fringing, a flux modulation with the wavelength that is characteristic 
of VLT/VISIR spectra and that has been found to be constant over periods of several months 
( Ivan Boekel et al.ll2009l ). This procedure may not be very accurate if the signal-to- noise ratio 
of the standard is very low, if the airmass of the standard is different from the airmass of the 
target, or if the photospheric lines of the standard are not well subtracted, because they are 
blended with telluric lines. To verify that none of these sources of error affect our results, 
we derive our calibrated spectra using an alternative approach. Specifically, we remove the 
fringing and the telluric lines using a fringing model and a model of the atmosphere and, 
then, we used the standard stars only for the flux calibratio n. A more detailed discussion of 
this alternative method is reported in lPascucci et al.l ( 1201 ll ). 



The two methods for deriving the calibrated spectra give the same results for all the 
stars except for T CrA. In this case, the first procedure did not work very well because the 
airmass of the star and of the standard are different and because the equivalent width of the 
[Ne II] line is very low, so systematic errors associated with flux calibration turn out to be 
larger than statistical errors. On the basis of these considerations, for T CrA we applied the 
second procedure for the flux calibration. 

Spitzer /IRS data for stars observed i n th e c2d survey were reduced usin g the c2d data 
reduction pipeline (see lLahuis et al.l 120071 and lLahuis &: Kessler-Silaccil 120061 for further de- 
tails). Data for other stars were downloaded from the Spitzer Heritage Archive in the basic 
calibrated data (BCD) format processed with the S18.18.0 data pipeline. These data were an- 
alyzed using custom IDL code (including SSC software IRSCLEAN MAS K) and the SMART 
(v. 8. 2.1.) program deve loped by the IRS Team at Cornell University (jHigdon et al.l 12004 ; 
Lebouteiller et al.ll2010l ). To further improve the flux calibration, we reduced 4 spectra of 
the standard star £ Draconis using the same procedure as for the science targets and, then 
we recalibrated the observed flux (the median among the 4 spectra) by using the median 
of these spectra and a model spectrum, that h as been produced by MARCS code with the 
specific purpose of calibrating Spitzer spectra ( IDecin et al.l 120041 ) . 



3. Results 

We observed 32 YSOs. Our targets belong to the star forming regions Taurus, p Oph, 
Corona Australis, Lupus and the Gum nebula and to the nearby associations TW Hya, e 
Cha and (3 Pic. 
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Our results are summarized in Table [3j Fluxes, full widths at half maximum (FWHMs) 
and peak velocities (v pea k) of the emission lines are calculated by fitting the spectral lines 
with the sum of a gaussian and a constant continuum (see Fig. [2]). Errors (la) in the 
total flux and FWHM are derived by a best fit procedure. Systematic errors associated with 
photometric calibration will be discussed later in this section. Errors on line peak velocities 
are the sum, in quadrature, of the statistical errors associated with the fitting procedure, 
the error on the stellar velocity reported in Table HJ and an a dditional error of 2 km s -1 that 



accounts for errors in positioning the star within the slit (see iPascucci et al.ll201ll and Sect. 
12. ip . If in the observed wavelength range there are no significant emission features, we assume 
that the line has not been detected and we report the upper limit in Table |3j The upper 
limit depends on the assumed line width and was computed as F up = 5a con V FWHMdv, 
where a con is the standard deviation of the continuum emission, FWHM is the assumed full 
width at half maximum of the line and 5v is the width of a velocity bin (~ 3.3 km s _1 ). 
Upper limits reported in Table [3] are calculated for an assumed FWHM=20 km s _1 . Four 
spectra of stars, for which [Ne II] was not detected are shown in Fig. [3] overplotted with a 
gaussian with a total flux equal to the upper limit and a FWHM=20 km s _1 . We chose to 
plot and report the upper limit for an assumed FWHM=20 km s _1 , because it is the typical 
FWHM of the observed emission line produced by photoevaporative winds (see Sect. 14.21) . 

Even if standard stars were observed just after or before the science targets and were 
chosen as close as possible to the targets, a residual systematic error due to discrepancy be 



tween the airmass of the standard and of the target could affect our results. IPascucci et al. 



(120111 ) studied the systematic errors on VISIR observations of the [Ne II] emission by observ- 
ing TW Hya several times during the same night. In their observations both flux and line 
peak velocity depended on the airmass. Specifically, the flux decreased of up a factor 2 as 
the airm ass increased from 1 t o 1.4, while the peak velocity varied from -4 to -7 km s _1 . In 
contrast, IPascucci et al.l ( 1201 ll ) found that the observed FWHM of the line did not depend 



on the airmass. 

All the detected lines are shown in Fig. [2] with best-fit gaussians overlaid on the data. 
All lines are well fitted by a gaussian, except RU Lupi and IRS 43. Specifically, the emission 
of RU Lupi looks asymmetric, with the line peak blue-shifted with respect to the line center; 
while in the spectra of both IRS 43 N and IRS 43 S we observe two features, a stronger one, 
which peaks very close to the stellar radial velocity and is well fitted by a gaussian profile, 
and a weaker one, which is strongly blue-shifted (v pea k ~ —100 km s _1 ) and asymmetric. 
However, this strongly blue-shifted component is very weak and we are not sure if it is a real 
feature or a statistical fluctuation. 



Two of our targets (Sz 73 and T Cha) have been previously observed by 



Pascucci fc Sterzik 
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( 120091 ). using VISIR. To test our data reduction procedure and t o compare our observation s 
with the previous ones, we rereduced again the data obtained by lPascucci fc Sterzikl ( 120091 ). 
using our version of the pipeline. We obtained slightly higher fluxes in both cases; these are 
the values reported in Table El The new version of the pipeline performs a better subtraction 
of the sky that can influence the observed flux of very faint targets like T Cha and Sz 73. 
However, these discrepancies are not significant, so they do not influence the scientific results 
of either ours or previous analysis. 

The two datasets for T Cha and Sz 73 also offer an opportunity to investigate variability 
of the [Ne II] line. In the case of T Cha the line flux, shape and FWHM remained the same. 



On the co ntrary, in Sz 73 we 



ailed to detect the [Ne II] from the outflow as previously 



reported in lPascucci &: Sterzikl (120091 ) . The second epoch data show upper limits of a factor 



1.6 lower than the epoch one data (see Table [3}. Because emission from Sz 73 is due to a 
protostellar jet, different slit orientations, different seeing or different positions of the star 
within the slit could produce significant discrepancies in the measured line flux. We can rule 
out the first two hypotheses, because in the two observations the slit orientation on the plane 
of the sky was the same and the seeing was slightly better during the second observation 
when the line was not detected. Errors in centering the star within the slit or intrinsic 
variability could explain the discrepancy between the results of the two observations. 

In Table O we also report line luminosities (or their upper limits), calculated from the 
distances reported in Table [TJ We also consider the effect of the absorption from interstellar 
and circumstellar gas, which, at 12.81 pm, is not ne gligible for highly emb edded sources, i.e 
when the absorption in the J band is A j > 1 mag ( iFlaccomio et al.l 120091 ) . Specifically, we 
obtained Aj (see Table |3]) from the literature and, when it was not available, from measure - 
ment of optical extinction, using the conversion law Aj = 0.282 Ay ( jRieke fc Lebofskylll985l ). 
or from the absorpt ion in the X-ray band (Nh), using the conversion law Aj = 1.8 x 10 _22 Nh 
( TVuong et al.ll2003l ). The absorption Aj is larger than 1 mag in 13 YSOs in p Ophiuchi, 3 
sources in the Corona Australis and, only marginally (Aj = 1.4 mag) in IRAS 08267-3336. 
Chapman et al.l (120091 ) showed that for highly absorbed stars in p Oph and other star forming 
region an extinction law with Ry = 5.5 is appropriate, as it accounts for grain growth in the 
dense part of the cloud. Therefore, for source s with A, T > 1, we calcu l ate th e absorption at 
12 //m, using the relation A 12 .si = 0.16Aj from lWeingartner fc Draind ( 1200 ll ). resulting in an 
absorption Ai 2 .8i ranging from 0.22 to 4.16 mag. For all other sources, we use d the standard 
R v = 3.1 extinction law (A12.81 = 0.097Aj from lWeingartner &: Draind 1200 ll ). However, for 
these sources, the absorption has a negligible effect (A 12 .si < 0.1 mag), independent of the 
extinction law used for the calculation. 



In Table [3j we also report [Ne II] fluxes measured from spectra taken by Spitzer /IRS 
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in the high resolution mode (R~600). We used values either reported from the literature or 
measured from archival data. For meas uring line fluxes or the ir upper limits from Spitzer 
data, we used the same procedure as in iFlaccomio et al.l (120091 ) . Specifically, we calculated 
the flux by integrating the emission between 12.78 and 12.84 /im and subtracting the underly- 
ing continuum by fitting a first or a second order polynomial in two adjacent spectral intervals 
(12.72-12.78 and 12.84-12.91 /an). In two cases (SSTc2dJ162145.1 and RX J1615.3-3255) 
for which high resolution spectra were not available, we used the low resolution spectrum 
and different wavelength ranges to estimate the line flux (12.7-12.9 /xm) and the continuum 
(12.3-12.7 and 12.9-13.3 urn). Upper limits (3 a) have been calculated, as previously dis- 
cussed for VLT/VISIR data, from the rms on the continuum emission and assuming a fixed 
line width (0.023 /jm). 



4. Discussion 

We detected the [Ne II] line in 12 YSOs (3 Class I, 4 Class II, 4 transition disks and 1 pre- 
transition disks); 11 of those are new detections with VISIR. We measured the FWHM and 
the velocity shift of all the detected lines. In Sect. 14.11 we compare fluxes of the [Ne II] line 
measured with VISIR and Spitzer and we discuss how discrepancies between observations 
with the two instruments can depend on the spatial location of the [Ne II] emitting gas. 
In Sect. 14.21 we compare all our results with other ancillary data (disk inclination, X-ray 
luminosity and mass accretion rate) and irradiated disk models to better understand the 
physical mechanisms generating [Ne II] emission from gas within the inner disk. 



4.1. Origin of the [Ne II] emission 



The aim of this work is to understand the origin of the [Ne II] emission in YSOs during 
the different phases of their evolution. Different emission mechanisms have been discussed 
in the literature. The [Ne II] emission could be a tracer of shoc ks produced by the interac 



tion between a protostella r outflow and the interstellar medium (jHollenbach fc McKedll989 



Hollenbach fc Gortil 120091 ) . This mechanism has been suggested to be a t the origin of the 
emission observed from the triple system T Tau (Ivan Boekel et al.l 120091 ). The main obser- 
vational evidence supporting this hypothesis is the spatial location of the emission: in fact, 
the emission is extended (1.1 ) and located along the stellar outflow. Otherwise, the neon 
could be ionized by high energy stellar emission. In this scenario the [Ne II] emission at 12.81 
fim should be produced within the inner circumstellar disk, few AU from the central star. 
A handful of observations of transition disks support this second scenario by showing that 
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the emission is usually compact and at a velocity close to that of the radial velocity of the 



centr al star ( IHerczeg et al.ll2007l ; iNajita et al.ll2009l ; iPascucci fc Sterzikll2009l ; iPascucci et al. 



2011( 1. 



Information on the spatial distribution of the [Ne II] emission can be obtained by com- 
paring fluxes measured by Spitzer (spatial resolution ~4.7 ) and VISIR (spatial resolution 
~0.4 ). To illustrate this, we show in Figure H] the ratio between the absorbed flux observed 
using Spitzer and VISIR as function of the [Ne II] absorption-corrected luminosity measured 
from the VISIR spectra for YSOs of various infrared Class. A systematic discrepancy is evi- 
dent between the [Ne II] fluxes measured by Spitzer and VISIR for Class I YSOs (blue dots); 
namely, Spitzer fluxes range from 2 to more than 20 times larger than [Ne II] fluxes measured 
by VISIR, with the exception of one object, which might thus be compatible with the other 
Class I stars. This discrepancy strongly suggets that the [Ne II] emission observed from 
Spitzer in most Class I objects is spatially extended, i.e. it is not produced within the inner 
disk (r<20-40 AU), but in the outer regions of the circumstellar envelope. Alternatively, the 
emission may span a very broad range of velocities, such that the upper limits reported in 
Table [3] could be underestimated (up to a factor 5 if we consider a line marginally resolved 
from Spitzer /IRS at high resolution mode). However, both scenarios would suggest that 
neon in class I objects is most likely ionized by shocks produced by protostellar outflows. To 
better investigate the [Ne II] emission from shocks located several AUs fro m the star, a dif- 



ferent observational strategy is required. As previously done for T Tau by Ivan Boekel et al. 



( 120091 ). the position and the velocity structure of multiple knots composing a jet and emit- 
ting in the [Ne II] can be studied by performing several long slit observations with different 
slit orientations to cover a large area around the target. Furthermore, to perform a complete 
study of the physical structure of protostellar jets, the results of these observations can be 
compared with narrow band imaging and spectroastrometric observations in other forbidden 
lines, like the [S II] doublet at 6717 and 6731 A or the [O I] at 6300 A. 

On the other hand, in YSOs which are classified as transition or pre-transition disks (Tr 
and pTr in Tableland black dots in Fig. H]), the ratio between [Ne II] fluxes measured by 
the two instruments is consistently less than 2, with the exception of two stars for which the 
upper limits might be compatible with other transition and pre-transition disks. For the 3 
objects detected with both instruments, the ratio is consistent with 1 at the 3 a confidence 
level. We conclude that, in transition and pre-transition disks, the bulk of the [Ne II] emission 
is produced within ~20-40 AUs of the central star. In the next section, we will discuss if 
this emission is due to a static disk atmosphere, a photoevaporative wind, or a magnetically 
accelerated outflow. 



It is more complicated to interpret results for Class II YSOs, namely stars harbouring 
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an optically thick disk. We detected the [Ne II] emission in 4 Class II YSOs. Since the ratio 
between Spitzer and VISIR fluxes ranges from ~1 to ~10, the spatially extended emission 
appears to be dominant, but significant emission from the inner disk may also be present. 



4.2. [Ne II] emission from the inner region: static disk atmosphere, 
photoevaporative wind or irradiated protostellar jet 

As discussed in the previous section, the detection of spatially unresolved [Ne II] emission 
in most of the transition and pre-transition disks demonstrates that the [Ne II] line is pro- 
duced within the inner disk (less than 20-40 Alls from the star) suggesting that neon is ionized 
by the high energy emission from the central star. How ever, the emission may be produced 
by a static disk atmosphere, a photoevapora tive wind (|Alexanderl 120081 ; lErcolano fc Owen 



20101 ) . or a magnetically accelerated outflow ( jShang et al.ll2010l ). The width and the velocity 



shift of the [Ne II] line at 12.81 /im should differ in these three cases. The emission from a 
static disk atmosphere should arise at the same mean radial velocity as the central star, with 
a broadenin g proportional to t he disk inclination (e.g. the [Ne II] emission from GM Aur 
observed by iNajita et al.l 120091 ) ; whereas a photoevapo rative wind - unless the disk has a 
large (~30 AU) inner disk hole (jErcolano fc Owenll2010l ) - will produce blue-shifted emission 
with respect to the central star (~ 1 — 10 km s~ x ) and, again, a broadening proportional to 
the disk inclination. Blue-shifted emission, very likely originating from a photoevaporative 
wi nd, has previously b e en ob served in 3 transitional disks (T Cha, CS Cha and TW Hydrae) 
by lPascucci fc Sterzikl (120091 ) . Finally, the emission line from a magnetically driven outflow 
should be blue-shifted to a larger velocity (~ 100 km s" 1 ) and muc h broader than emission 
from a disk, as observed, e.g. in Sz 73 by lPascucci fc Sterzikl (120091 ) . 



In Fig. |5l we show the FWHMs of the detected lines as function of their blue-shifts with 
respect to the central star velocity. Colors represent the different infrared Classes, as in Fig. 
IU Our results rule out the possibility that the [Ne II] emission is produced by a static disk 
atmosphere, since all of the detected lines are blue-shifted with respect to the central star. 
In two of the stars of our sample (RU Lupi and IRS 60) the emission likely originates from 
a magnetically driven outflow, as demonstrated by the large blue-shift (v pe afe = 168±4 and 
Vpeafc = 55±6 km s -1 , respectively) and the large FWHM (107±9 and 49±12 km s _1 , respec- 
tively) of the observed emission lines. This hypothesis is confirmed by o bservational evidenc e 
that indicate RU Lupi is the source of a well-known protostel lar jet (e.g . iTakami et al.ll200ll ). 
and by the proximity of IRS 60 to a Herbig-Haro objects (jWu et al. 20021 ). Furthermore, 
the line profile of RU Lupi (top-right panel of Fig. [2]) appears slightly asymmetric, with 
the bluer part of the profile decreasing faster than the redder part. Similar asymmetric 
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and strongly blue-shifted profiles have been predicted by IShang et al.l (120 101 ). who modeled 
the [Ne II] emission produced by a magnetically accelerated wind irradiated by stellar high 
energy emission. Unfortunately, the signal-to-noise is too low to carry out a detailed com- 
parison between our spectra and these models. In the spectra of IRS 43 N and IRS 43 S, we 
also observe a secondary emission component that is strongly blue-shifted and asymmetric. 
Ho wever, although in each case the profile of this component resembles the profiles predicted 
by lShang et al.l (120101 ). these highly blue-shifted and asymmetric components from both IRS 
43 N and IRS 43 S are too weak to establish if the resemblance is significant. 

The emission observed from the other 10 YSOs is blue-shifted by less than 18 km s -1 . 
For 5 of our targets, specifically transition and pre-transition disks, a disk origin is confirmed 
by the correlation between line FWHM and the disk inclinatio n. The correlation is show n 
in the left panel of Fig. El where we also plot two points from iPascucci fc Sterzik! ( 120091 ). 
In this panel, we also report the relation between FWHM and disk in c linati on derived from 
irradiated primordial and transition disk models by lErcolano fc Owen! ( 120101 ) . Our measure- 
ments agree very well with the models, except in the case of RU Lupi and T Cha (the points 
on the top-left and top-right of the plot). In the first case, the emission is likely not due 
to a photoevaporative wind, but rather to a mag netically accelerated outflow, as discussed 
above, so the models of lErcolano &: Owen! (120101 ) are not applicable. In the case of T Cha 
the discrepancy between models and obs ervational data may be related to the difference 
between the mass of T Cha (1.3 M Q from ISchisano et al.l 120091 ) and the mass of the central 
star assumed in the models (0.7 M Q ). Discrepancies between the models and observations 
are more evident in the right panel of Fig. [HI where we plot the line blue-shift as a function 
of disk inclination. Specifically, the obser ved blue-shi f ts are higher than the predicted ones. 
Higher blue-shifts have been predicted by lAlexanderl (120081 ). who considered a warmer disk 
atmosphere (T~10,000 K) irradiated by stellar EUV radiation. 

For the other 5 YSOs with an observed blue-shift of less than 18 km s -1 (IRS 43 N, 
IRS 43 S, IRS 45, SSTc2dJ162145. 13-234231.6 and T CrA), the origin of the emission is less 
clear, because we do not know the disk inclination. Furthermore, four of these objects are in 
the earliest stage of the disk evolution (3 Class I IRS 43 N, IRS 43 S and SSTc2dJ162145.13- 
234231.6 and a Class II IRS 45) and the line FWHM is larger than observed in the transition 
disks, so the [Ne II] line at 12.81 /im could be produced by a jet seen with a large angle 
respect to the line of sight. We also note that in these objects errors on the blue-shift could be 
underestimated, because no radial velocities measurements are available in the literature, so 
we based our radial velocity estimate on the radial velocity of the parent cloud (p Ophiuchi) 
or other nearby stars (see Tabled]). Observations aimed to measure accurate radial velocities 
and disk inclinations for these stars can help to distinguish between jet and disk origins. 
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In the left panel of Fig. [TJ we report the [Ne III lumin osity as a function of the X-ray 
luminosity, with the predictions of lErcolano fe Owenl (120101 ) overplotted with a dashed line. 
[Ne II] emission and X-ray luminositi es do not appear correlated, as would be exp ected 
from the irradiated disk models (e.g. Meiierink et al.l 12008k lErcolano fc Owenl |2010| ) and 
as has been observed by Giidel et al. ( 201ol ). However, the lack of a correlation between 
[Ne II] emission and X-ray luminosities could be due to the non-homogeneity of the observed 
sample. Specifically, our sample includes very different disks around stars of different mass 
that are going through different phases of the star formation process, while irradiated-disk 
models that predict a correlation focused their analysis on a single prototypical disk. The 
relation between disk properties and predicted [Ne II] luminosity has been investigated by 



Schisano et al.l (120101 ). who found that varying disk properties can cause luminosity variations 



up to one order of magnitude. Several up per limits in Fig. [7] (left p anel) lie below the [Ne II] 
luminosity predicted from the model of lErcolano fc Owenl (120101 ) . The lack of detectable 
[Ne II] emission in Class I and Class II sources could depend on the effects of absorption on 
X-ray and EUV emission. Specifically, if EUV and X-ray emission are strongly attenuated by 
the optically thick inner disk layers, the amount of ionized neon could be less than predicted 
by the models. In transition disks - like Hen 3-600, for which the upper limit on the [Ne II] 
luminosity is one order of magnitude below the expected value - the discrepancy between 
model and observations may be attributed to the opposite reason; namely, some of these disks 
have lost a large part of their circumstellar gas, so we do not see strong emission lines in the 
mid-infrared. In the case of Hen 3-600, a non-detecti on of [Ne II] emission is consistent both 
with the lack of detection of molecular CO emission (I Andrews et al.ll2010l ) and its low mass 
accretion rate (se e Table [H). In contrast, higher mass accretion rate s and strong molecular 
CO emission (see Kastner et al. 1997 . 20081 . 2010 ; Curran et al. 2011 ) characterize TW Hya, 
V4046 Sgr and MP Mus where a strong [Ne II] emission was detected. Another possibility 
is that neon is ionized by EUV or soft X-ray photons mainly produced by accretion shocks, 
hence the line emission a 12.81 /iin strongly decreases at low mass accretion rates. 

The right panel of Fig. [7J shows a weak correlation between [Ne II] luminosity and mass 
accretion rate, but this correlation is strongly influenced by just three objects (RU Lupi, IRS 
43 N and IRS 43 S). As discussed above, the [Ne II] emission in RU Lupi is produced by a 
protostellar jet, wh ile the origin of the em issio n from IRS 43 N an d IRS 43 S is uncertain. 
So, as suggested by lFlaccomio et al.l (120091 ) and lGiidel et al.l (120101 ) . the correlation between 
[Ne II] luminosity and mass accretion rate may be driven by the presence of strong emission 
from jet sources, which are typically strong accretors. 
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5. Summary and Conclusions 

We observed 32 YSOs belonging to different infrared Classes using the high resolution 
(R~30000) mid-infrared spectrograph VISIR at the VLT, with the aim of studying the origin 
of the [Ne II] emission line at 12.81 fim. We detected the line in 12 YSOs, thus tripling the 
number of detections of 12.81 fim [Ne II] emission in YSOs at high spectral resolution. 
We also collected from the literature published Spitzer /IRS data for 14 of our targets and 
analyzed archival data for another 14 stars. By comparing VLT/ VISIR and Spitzer /IRS 
data with ancillary data retrieved from the literature, we obtained the following results: 

• Comparison between VISIR (slit width 0.4 ) and IRS (slit width~ 4.7 ) fluxes demon- 
strates that for Class I YSOs the [Ne II] emission is mainly produced by gas located at 
a distance of more than 20-40 AU (i.e., outside of the VISIR slit but within the IRS 
slit) from the central star and, therefore, arises from the extended envelopes. At these 
distances from the star, neon is most likely ionized by energetic shocks produced from 
protostellar outflows. 

• The same comparison (VISIR vs. IRS) shows that the [Ne II] emission from transition 
and pre-transition disks is produced by gas located within a few AU from the central 
stars, while in Class II objects we observe emission both from the inner region and 
from the extended envelope. 

• Observed blue-shifts and FWHMs of the [Ne II] line are consistent with emission in 
either a photoevaporative wind or the launching region of magnetically driven outflow. 
We associated the emission with a photoevaporative wind in 5 pre-main sequence 
stars with transition and/or pre-transition disks (blue-shifts between 2 and 18 km s" 1 , 
FWHM between 15 and 45 km s -1 ), where we observe a correlation between line 
FWHM and disk inclination. Otherwise, we associated the emission to an outflow in 
two objects (blue-shifts of 55 and 168 km s _1 , FWHM of 49 and 107 km s _1 ), one of 
which, RU Lupi, is a prototypical target for studies of protostellar jets. The origin of 
the emission in the remaining 5 objects is more uncertain. 

• We compared the dependence of line FWH M and blue-shift on disk i nclinations with 



the predictions of irradiated disk models of lErcolano Sz Owed (120101 ) . The observed 
relation between line FWHM and disk inclination is consistent with the predictions of 
models, except for a single star (T Cha) whose disk is seen at high inclination angle. 
Line blue-shifts are larger than predicted from the models. 

Blue-shifted emission and a correlation of the line width with disk inclination have 
been also predicted in the case of photoevaporative wind triggered by EUV radiation 
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([Alexander! 120081 ). 



• We detected [Ne II] emission with a blue-shifts smaller than ~ 20 km s _1 from 3 Class I 
objects and 2 Class II objects. Although these projected velocities appears small for an 
outflow origin, and would rather point toward a disk wind, we have no constraints on 
the inclinations of these systems and thus on the true velocity of the emitting material. 

Our observations demonstrate that the [Ne II] emission in YSOs arises either from 
shocks formed in protostellar jets or from a photoevaporative wind in the inner disk, and 
high spatial and spectral resolutions observations are necessary to distinguish between the 
two mechanisms. Specifically, the former mechanism prevails in early stages of the star 
formation process when powerful outflows are generated, while the latter prevails in the 
latest stages of disk evolution. However, to understand if the two processes may coexist 
and photoevaporative winds are already active at early stages, a better knowledge of disk 
properties (i.e. disk inclination, radial velocities and disk structure) is required. The new 
instruments operating in mid-infrared and submillimeter band (VLT/CRIRES, Herschel, 
ALMA) will allow us to better study these systems and interpret our results. 

For what concerns the more evolved transition and pre-transition disks, the presence 
of emission blue-shifted of ~ 2 — 10 km s _1 and the correlation between line width and 
disk inclination are strong evidences for the photoevaporation from the disk. However, 
observations in the [Ne II] line are not sufficient to understand the full physical scenario, 
specifically, it is not clear if EUV emission produced by accretion shocks or X-rays produced 
either by accretion shocks or coronae is the main driver of photoevaporation and disk heating. 
High spect ral resolution observa tions of other forbidden forbidden emission lines (e.g. [O I] 



at 6300 A, IPascucci et all l201ll ) emitted in the inner disk region may help to address this 
issue. 

Finally, we detected the [Ne II] line in only 13 YSOs out of the 32 observed. This is 
not surprising given that upper limits are in many cases of the same order of magnitude of 
the expected luminosity (see left panel of Fig. [7]). Therefore, our study is biased toward 
stronger and closer sources. A higher sensitivity instrument is required to observe a larger 
and more complete sample and thereby to fully explore the space of physical parameters of 
the star-disk system. The upgraded VLT/VISIR expected in the next summer will be a first 
step in this direction. 
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iKenvon et"al] (1199411 ; (2) iFurlan et ail J20091) ; (3) lHartmann et all fa98Tb ; (41 lEspaillat et all ll2007al l; (5) IPetterssonl (120081 1 ; (61 iGiidel et all d2010h ; (71 iKastner et all 
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a Full stellar name SSTc2dJ162145.13-234231.6. 

b Mean radial velocity of ten memebers of the p Oph star forming region measured bv lDoppmann et al. I i2003l) . 

c The heliocentric velocity has been derived from the velocity with respect to the local standard of rest using the IRAF task rvcorrect. 

d Assumed equal to the radial velocity of IRS 47. 

°Measured by submillimeter interferometric observatins of the disk. 

f Derived from the fit of the SED. 

s Measured by scattered-light imaging of the disk. 

h Assumed equal to the inclination of the stellar rotation axes. 

'Unabsorbed X-ray luminosity in the range 0.3-10 KeV. 
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Table 2. Observation log. 
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Table 2 — Continued 
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Table 3. Results. 
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10 — 14 erg cm~ 2 s~ 1 km km s — 1 mag 10 28 erg s^ 1 10~ 14 erg cm~ 2 s~ 1 



1 


V807 Tau 


<1.1 


- 


- 


0.2 


< 2.6 


- 


1 - 


2 


LkCa 15 


<0.5 


- 


- 


0.3 


< 1.2 


0.28± 0.02 


1 11 


3 


IRAS 08267-3336 


<0.5 


- 


- 


1.4 


<13.8 


1.40± 0.08 


4 11 


1 


Hen 3-600 A 


<0.2 


- 


- 


0.2 


<0.03 


< 0.5 


2 11 


5 


WX Cha 


<0.2 


- 


- 


0.6 


<0.03 


0.57± 0.06 


1 11 


6 


XX Cha 


<0.2 


- 


- 


0.3 


<0.03 


0.48± 0.02 


1 11 


7 


T Cha 


3.4±0.3 


44.9 ±3.2 


-10.5±2.7 


0.8 


5.1±0.5 


3.20± 0.21 


5 12 


7 


T Cha a 


3.1±0.2 


39.4±1.9 


-4.7±2.5 


0.8 


4.7±0.3 


3.20± 0.21 


5 12 


8 


MP Mus 


l.lzfcO.l 


15.9 ± 1.4 


-4.4±2.1 


0.1 


1.4±0.2 


< 2.0 


2 13 


9 


IM Lupi 


<0.5 


- 


- 


0.2 


< 1.4 


1.07± 0.05 


4 11 


10 


SZ 73 


<0.6 


- 


- 


0.8 


< 1.8 


1.60± 0.24 


6 12 


10 


Sz 73 a 


1.9±0.2 


72.8±6.4 


-97±4 


0.8 


5.4±0.6 


1.60± 0.24 


6 12 


11 


RU Lupi 


2.4±0.2(<0.3) b 


106.7 ±8.8 


-168.2±4.0 


0.0 


6.4±0.7 


2.60± 0.83 


2 4 


12 


RX J 16 15. 3-3255 


1.4±0.2 


20.5 ± 2.7 


-7.5±2.8 


0.3 


3.8±0.7 


2.76± 0.46 


7 11 


13 


SSTc2dJ162145.1 


5.8±0.5 


44.1 ± 2.6 


-18.0±2.7 


2.8 


15.0±1.2 


11.83± 0.46 


8 11 


14 


DoAr 25/GY92 17 


<0.3 






0.7 


< 0.5 


0.50± 0.10 


3 3 


15 


WL 10/GY92 211 


<0.1 






4.5 


< 0.4 


1.60± 0.30 


3 3 


16 


Elias 29/GY92 214 


<0.5 






11.4 


< 4.8 


< 8.5 


3 3 


17 


SR 21 


0.5±0.1 


15.1 ± 1.2 


-8.3±2.7 


2.3 


1.3±0.1 


< 3.0 


1 4 


18 


WL 19 


<0.8 






16.3 


<15.2 


< 0.8 


3 3 


19 


WL 20/GY92 240 


<0.2 






4.1 


< 0.7 


6.28± 0.25 


3 3 


20 


IRS 43/GY92 265 S 


3.7±0.3 


37.5 ± 2.1 


-6.2±2.7 


8.2 


21.2±1.5 


45.30± 1.70 


3 3 


21 


IRS 43/GY92 265 N 


1.6±0.2 


35.7 ±4.3 


-7.6±3.1 


8.2 


9.1±1.4 


45.30± 1.70 


3 3 


22 


IRS 44/GY92 269 


<0.7 






12.8 


< 8.2 


8.00± 2.40 


3 3 


23 


IRS 45/GY92 273 


0.6±0.1 


56.4 ±6.4 


-12.4±3.7 


6.6 


2.9±0.4 


2.19± 0.41 


3 3 


24 


IRS 47/GY92 279 


<0.2 






7.4 


< 1.0 


2.01± 0.51 


3 3 


25 


IRS 60/WSB 71a 


0.9±0.3(<0.2) b 


49.4 ± 12.0 


-54.9±5.7 


2.5 


2.3±0.8 


1.91± 0.09 


1 11 


26 


DoAr 44 


<0.3 






1.2 


< 0.5 




8 - 


27 


V4046 Sgr 


6.6±0.2 


22.5 ±0.5 


-10.5±2.0 


0.0 


4.2±0.1 


7.68± 0.22 


2 11 


28 


RX J1842.9-3532 


<0.2 






0.3 


< 0.5 


0.43± 0.13 


9 13 


29 


CrA IRS5 A 


<0.6 






7.0 


< 3.3 


20.32± 0.65 


10 11 


30 


CrA IRS5 B 


<0.6 






7.0 


< 3.1 


20.32± 0.65 


10 11 


31 


R CrA 7B 


<0.5 






26.0 


<42.9 


10.51± 0.42 


10 11 


32 


T CrA 


0.7±0.2 


19.1 ±4.7 


-2.0±3.0 


0.7 


1.6±0.5 


8.88± 1.30 


10 11 



Note. — For the Ne II detections, fluxes, FWHMs and the line peak velocity (v pea fe) are calculated by fitting the spectrum with a gaussin 
profile for the line and a first order polynomial for the continuum. The line peak velocity is measured with respect to a reference system 
comoving with the star. Errors (ltr) are derived from the best-fit procedure. The Ne II luminosity are corrected for the extinction as discussed 
in the text. Th e last column repor ts t he references for the Aj and the Ne II flux as meas ured from Spitzer low and medium resolutio n 
spe ctoscopy: (1) Furlan et al l [|2009h; (2l|Curran et all [|201ll); (31 iFlaccomio et al.l J2009f); (4) iGiidel et afllboioll; (51 ISchisano et all 20091') ; 
(61 iHuehes et al [|l994ll; (71 Merfn et alj fcoicft; (8) lEyans et all J2009h ; (91 ICarpenter et alj feOOSll ; (10) iForbrich fc Preibischl ||2007I 1; Till 
this work; (12) [ 



Lahuis et al 



200 



3); (13) 



Pascucci et alj (2007) 



a Results obtained from a new reduction and calibration of the observations carried out by Pascucci & Stcrzik ( 2009) 



'Upper limit for a disk emission centered at km . 
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Fig. 1. — Examples of reduction of VISIR spectral image data. The left panels show weight 
maps calculated by the pipeline, while the central and the right panels show the final 2D 
spectrum obtained by the pipeline and by an alternative data reduction using our own 
weight map, respectively. From top to bottom we show the frames for a very bright star (the 
standard HD 151680), a very faint star (WL 10) and a binary (CrA IRS 5), respectively. 
The red and the green lines overplotted on the left panels are the profiles along the spatial 
direction of the weight maps calculated from the pipeline (red line) and defined by us (green 
line). The two profiles are normalized to the same maximum value. On the x and y axes, 
we report the number of pixel and the intensity of the profile of the weight, respectively. For 

"f" Vi o KmQTMr Tiro o Vi r\nr r^nlir "f" Vi o -XTCTtui rrln-f - TnoT~\ "fv \ v nonnn^ovir r~> r\ Tin ri ati onf 
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Fig. 2. — Flux calibrated spectra of the 12 sources, for which we detected the 12.81 /im 
[Ne II] emission line. The gaussian best fit of the data is shown by the red curves. The 
spectra are plotted as function of radial velocity in the stellocentric reference system. 
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Fig. 3. — Flux calibrated spectra of the 4 sources for which we did not detect the 12.81 /zm 
[Ne II] emission line. Overplotted red lines show gaussian profiles of a line with a total flux 
equal to the upper limit and a FWHM=20km s™ 1 . Radial velocities (x axes) are with respect 
to the observed star as in Fig. [2j 
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Fig. 4. — Ratio between flux of the [Ne II] line as measured from Spitzer /IRS and 
VLT/VISIR spectroscopic data as function of the [Ne II] luminosity measured from the 
VISIR data. Blue and red dots indicate Class I and Class II YSOs, respectively, while tran- 
sition and pre-transitional disks are both indicated with black dots. For the two binary 
systems (IRS 43 and CrA IRS5) not resolved by Spitzer, but resolved by VISIR, we assumed 
that the [Ne II] emission observed by Spitzer was equally shared by the two components. 
Numbers are the IDs of the stars reported in Table [3j 
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Fig. 5. — FWHM of the [Ne II] emission line as function of the blue-shift with respect to 
the stellar velocity. Colors of the symbols are as in Fig. |H Numbers are the IDs of the stars 
reported in Table [3] 
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Fig. 6. — FWHM (left panel) and blue-shift (right panel) of the [Ne II ] emiss ion as functions 
of the disk inclination. Empty symbols are from iPascucci fc Sterzikl (120091 ) . filled symbols 
from this work. Colors of the sy mbols are as in Fig . Pfl In both panels, we overplot theoretical 
predictions from the models of lErcolano fc Owen! (120101 ). Specifically, red, green, blue and 
cyan continuous lines were computed for transition disks irradiated by a star with X-ray 
luminosity Lx = 2 x 10 30 erg s _1 and with an inner hole of 8.3, 14.2, 21.1 and 30.5 AUs, 
respectively. Magenta, orange and black dashed lines were computed for a primordial disk 
model without inner hole irradiated by a star with an X-ray luminosity of 2 x 10 28 , 2 x 10 29 
and 2 x 10 30 erg s^ 1 , respectively. Numbers are the IDs of the stars reported in Table [3] 
(object 11 is RU Lupi). 
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Fig. 7. — Luminosity in the [Ne II] line as function of the X-ray luminosity in the 0.3-10 
keV range (left panel) and mass accretion rate (right panel). Colors of the symbols are as 
in Fig. HI while the empty symbol shows the upper limit for RU Lupi for a second emission 
component located at zero velocity. The X-ray luminosity of the binary system IRS 43 has 
been assumed to be equally shared by the two components. The dashed line on the left panels 
represents the relation between [Ne III arid X-ray luminosity estimated from the irradiated 
disk models of lErcolano &: Owen! (120101 ) . Numbers are the IDs of the stars reported in Table 

El 



